Mon. Not. R. Astron. Soc. OOO,[TH6](2011) Printed 8 May 2012 (MN WT^ style file v2.2) 



Relativistic iron K X-ray Reverberation in NGC 4151 



(N 
O 



O 
(N 

> 



X 



A. Zoghbi^'^*, A. C. Fabian^, C. S. Reynolds^ and E. M. Cackett^'^ 

^Department of Astronomy, University of Maryland, College Park, MD 20742, USA 
'^Institute of Astronomy, Madingley Road, Cambridge CBS OH A 

^Department of Physics and Astronomy, Wayne State University, Detroit, MI 4-8201, USA 



ABSTRACT 

Recent X-ray observations have enabled the study of reverberation delays in AGN for 
the first time. All the detections so far are in sources with a strong soft excess, and the 
measured delay is between the hard (1-3 keV) direct continuum and the soft excess 
(0.5-1 keV), interpreted as the reflection continuum smeared by relativistic effects. 
There is however an inherent ambiguity in identifying and studying the details of the 
lines in the soft excess. Here we report the first detection of reverberation in the iron 
K band in any AGN. Using XMM-Newton observations of NGC 4151, we find delays 
of order 2000 s on time-scales of 10^ s between the 5 — 6 keV band and 2 — 3 and 7 — 8 
keV bands, with a broad lag profile resembling a relativistically-broadened iron line. 
The peak of the lag spectra shifts to lower energies at higher frequencies, consistent 
with the red wing of the line being emitted at smaller radii, as expected from reflection 
off the inner accretion disk. This is a first detection of a broad iron line using timing 
studies. 

Key words: X-rays: galaxies - galaxies: individual: NGC 4151 - galaxies: active - 
galaxies: Seyfert - galaxies: nuclei. 



• ■ 1 INTRODUCTION 



X-ray observations of active galactic nuclei (AGN) have 
shown that most of the radiation is emitted from very close 
to the central super-massive black hole. Broadening of X- 
ray spectral features caused by relativisti c effects has been 
the main tool in studying these systems (Milleij [20071 for a 
review). Recent work has moved beyond spectral modelling 
to detecting time delays between the different spectral com- 
ponents. Short time delays have first be en seen in the Nar- 
row Line Seyfert 1 galaxy 1H0707-495 (|Fabian et al.ll2009l : 
IZoghbi et ailboiol ). Thanks to its strong soft excess and high 
iron abundance, small delays of a few tens of seconds are de- 
tected between the 1-4 and 0.5-1 keV bands where the soft 
band lags the hard band (hence soft lag) . The natural expla- 
nation for the lag is reprocessing by the reflecting accretion 
disc in the vicinity of the black hole . Similar delays have also 
been seen in several other objects jEmmanoulopoulos et al.l 
I2OIII : Ide Marco etHlbOUl : IZoghbi fc Fabianll201ll ). A com- 
mon feature among all these sources is a strong soft excess in 
the spectrum below 1 keV. This is interpreted, in line with 
the lags, as a relativistically-blurred reflection continuum, 
that responds to variations in a primary illuminating hard 
power-law. 

The interpretation of the soft lag (soft excess lagging 
the hard continuum) as a light crossing efl^ect close to the 
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black hole was questioned bv lMiller et al.l l|2010l ). who sug- 
gested that such small delays could be an artefact of a much 
larger reverberating system of clouds (~1000 instead of ~ 
a few gravitational radii Vg ). Although this interpretation 
fails when a complete set of observa bles (variability pow er, 
flux and lag spectra) are considered (jZoghbi et al.|[201ll ). it 
stan ds on the ambiguity in interpreting the s oft excess itself 
(e.g. ICrummv et"al] I2OO6I : [DoneliaiJ[20iJ). The smooth- 
ness of its spectrum makes identifying the relativistically- 
broadened emission lines therein difficult unless the object 
has some extreme properties (e.g. high iron abundance and 
strong Fe-L in 1H0707-495). Therefore, a key challenge for 
any model is its ability to explain the reverberation in the 
iron K region (6-7 keV) where such an ambiguity is not 
an issue. The vital task is to find and observe an object 
that has reverberation signatures in this band. The impor- 
tance of reverberation studies in the K band have been dis- 
cussed previously, along with theoretical calculations of the 
expected behaviour of the broad iron line with t ime (jStellal 
I1990I : iMatt fc Peroialll992l : [Reynolds et"allll999l ). However, 
it remained difficult to detect these reverberation delays in 
real data. 

The count rate at the bands of interest is a key pa- 
rameter in reverberation studies. All previous lag detections 
have been in the soft band (< 1 keV) in objects with a 
strong soft excess. Those objects however, show steep spec- 
tra where the flux drops significantly in the iron-K band. 
Therefore, a different class of objects need to be explored if 
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Figure 1. Unfolded spectra (in EF{E) form) of NGC 4151 from 
the XMM-Newton PN detector. The labels include: Obs. number 
(Obs. ID : total exposure time). 
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Figure 2. Left: Lag spectrum at frequencies smaller than 5 X 10~* 
Hz for NGC 4151 using all the eight observations. The delays are 
measured with respect to the total 2-10 keV band (i.e. the lags of 
all the points are measured with respect to roughly their average). 
Right: Same as Left but only for the lowest frequencies. 



Fe-K reverberation is to be studied. In this work we show 
that NGC 4151, which is ~ 300 times brighter than 1H0707- 
495 at 6.4 keV, shows inter-band delays that indicate that a 
significant fraction of the flux in this band originates in the 
vicinity of the black hole. 

NGC 4151 is a well-studied Seyfert galaxy {z = 0.0033) 
that harbours one of the brightest AGN in the X-ray band. It 
has been observed and studied by all major X-ray satellites. 
The X-ray to 7-ray spectrum is dominated by a Comptonisa- 
tion component that breaks at ~ 100 keV. There is a strong 
narrow iron line at 6.4 keV and a r e flection hump peak- 
ing a t 30 keV l|Zdziarski et all Il996l : ISchurch fc Warwickl 
|2002| ). These components are absorbed below 4 keV by a 
high column {Nh ^ lO^^cm"^) of gas in the line of sight 
l|Weaver et al.lllS)94l ) . The spectrum below 2 keV is very com- 
plex. High resolution Chandra observations revealed that 
it is spatially-resolved and it is associated with a highly 
ionise d plasma c oincident with the opti cal narrow-line re- 
gion l|Ogle et al.|[2000 ; Wang et al.ll201ll ). This component 
does not vary much and will not be discussed here. 

NGC 4151 has been observed with XMM-Newton sev- 
eral times in short observations constrained by its visibil- 
ity. The unfolded spectra are shown in Fig. [1] for the eight 
longest observations (~ 19 — 62 ks) from a total of 13 XMM- 
Newton. There is clearly significant spectral variability. The 
difference between the spectra can be as simple as a uniform 
apparent increase in flux, or it could be due to changes in 
the absorber on time scales of years (e.g. difference between 
observations 1 and 7 in Fig. [1} . 

In this work, fight curves from eight archival XMM ob- 
servations are used (Obs. IDs and length of observations 
are shown in Fig. [T|. ODF data were reduced using SAS 
V. 11.0.0. Light curves are produced in equal energy steps 
above 2 keV with 1024 s time bins. Instrumental corrections 
were applied using the task epiclccorr. The lag was then 
estima ted at the observed times following iRybicki fc Presa 
l|l992h . We determine the uncertainty in the lags via Monte 
Carlo simulations. We shift each point in the hght curves 
by drawing a random number from a Gaussian distribution 
with mean equal to the count rate and standard deviation 
equal to the uncertainty in the count rate. This is performed 
100 times, and the 1-sigma uncertainty in the lag is taken 
from the resulting distribution. The errors are very simi- 



lar t o those estimated in the standard way l|Nowak et al.l 
I1999I) . Lag-energy plots are produced by calculating the lag 
at the frequency of interest between the energy band and 
the fight curve in the whole 2-10 keV band excluding the 
current band. This maximises the signal to noise and es- 
sentially measures the delays relative to the average (e.g. 
IZoghbi et aLll201ll ). 



2 LAG MEASUREMENT IN NGC 4151 

The mass of NGC 4151 is determined accurately with 
reverberation ma.p ping to be AI = 4.5 ± 0.5 x 10^ 
M0 (|Bentz et al.ll2006l ). If the time scales seen in 1H0707- 
495 and RE J1034-I-396 are scaled with mass, reverberation 
delays are expected at ~ 3 x 10~^ Hz. Fig. [2]-left shows the 
lag spectrum for all the observations measured at frequen- 
cies < 5 X 10~"^ Hz , calculated in a similar way to those in 
IZoghbi et all (|201ll ). 

In the right panel of Fig. (2] we show the lag spectra for 
the lowest frequencies (derived from a single frequency bin 
estimate) as they are estimated from the longest observa- 
tion (i.e. the 2nd shown in Fig. [l}. This plot clearly shows 
that the ~ 6 keV band lags the > 8 and < 4 keV bands 
by ~ 2000 seconds on time scales of ~ 10^ seconds. Since 
this observation happens to map a low flux source state this 
could be a possible indication of a relation between the lag 
and the flux state. 

Calculating the lags for the low and high flux intervals 
separately shows that the lags are prominent only in the low 
flux observations. To avoid the obvious effects of changes 
due to the absorber, low flux observations include observa- 
tion 1,2 and 3 (as defined in Fig. [TJ. The high flux observa- 
tions are 4, 5 and 6. Fig.|3]-top shows the lag-energy plot for 
the two flux intervals calculated between (0.5 — 5) x 10~* 
Hz. The lower limit is dictated by the length of the high 
flux observations. The conclusions are not very sensitive to 
the upper limits apart from a scaling factor for the lag and 
slightly larger errors if smaller band is considered. Again, at 
low fluxes, the ~ 5 keV band lags both the 2 — 3 and 8 — 9 
keV bands, while the lag spectrum is significantly different 
at high fluxes. The plot clearly shows an interesting struc- 
ture at the iron K energies that resembles a broad iron line. 
It should be noted here that the lags are not very promi- 
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Figure 4. Lag vs frequency plot for the low flux interval. Hard 
lags are plotted, so a positive delay indicate hard bands lagging 
softer bands. Left: The lag between 2—3 and 4—5 keV bands. Right: 
The lag between 2—3 and 5—6 keV. The energy and frequency 
dependence of the lag is apparent. 



Figure 3. Top: Lag spectrum at frequencies (0.5 — 5) X 10~^ Hz 
for NGC 4151 for the low and high flux states separately. The 
plots are similar to Fig. [2] Lags for the low flux have been shifted 
up by 1500 seconds for visual clarity. Bottom: Reflection fraction 
calculated as the ratio of the relativistically blurred component 
and the power-law component taken from the best fit to observa- 
tions 1 and 4 representing low and high flux intervals respectively 
(see SeclSj. Vertical lines mark the peaks in both the lag and flux 
spectra (at 4.9 and 6.4 keV for the lag spectrum and reflection 
fraction respectively) . 



nent when all observations are included because the high- 
flux signal dominates when combining the data. It should 
be noted also that the peak in the lag spectrum measured 
at (0.5 — 5) X 10~* Hz is at lower energies compared to that 
measured at lower frequencies (Fig. [3]-top vs. [2]), an effect 
which is discussed in Sec|3] 

The frequency-dependent lag is shown in Fig. |4] for the 
low flux observations. It shows how the lag depends on en- 
ergy as well as frequency. The lag peaks at lower frequen- 
cies when higher energies are considered (Fig. |4]-right), and 
the peak is at higher frequencies when a lower energy band 
is considered (Fig. [IJleft). Due to the length of the obser- 
vations, the lowest frequency point contains one frequency 
value (from the longest observation). This can be problem- 
atic because of the small number of points included in the 
averaging process and only longer future observations will 
provide any improvements. There are however several rea- 
sons that strengthen our measureme nts of an ene r gy an d 
frequency dependence. As shown in IZoghbi et al.l l(20ld ). 
Monte Carlo simulations show that the lowest frequency 
point always shows a bias toward zero and the bias is pro- 
portional to the absolute value of the lag (when measured 
at the same frequency). Therefore, because the difference 
in value of the lowest frequency point between the left and 
right of Fig. |4]is significant, their values might be biased but 
their difference has to be intrinsic. This is further supported 
by a clearer difference at the two higher frequency points, 
where the bias is irrelevant. 



3 LAG INTERPRETATION. 

To understand these lag spectra, we turn to the flux spec- 
tra. It was realised from early observations that the source 
is seen through a high column density absorber. How- 
ever, a model consisting of a highly absorbed power-law is 
not sufficient to fit the spectral curvature in the 2-5 keV 
band (jWeaver et al.lll994l ). Instead, it can be reproduced by 
adding a partially ionised (warm) absorber, another neutral 
absorber, or a combinat ion where some abs o rbers only par- 
tially cover the source (Weaver et al." 1994; Warwic k et al.l 
[1995; S church fc Wa rwick 2002; Puccotti et al. 2007). How- 
ever, in a systematic analysis of the Fe K band for a sample 
of AGN, iNandra et all (|2007l ) noted that there is red wine 
to the line that cannot be accounted for by the absorbers. 
A similar conclusion was reached earlier by lYagoob et al.l 
(1995>). 

The lags shown in the top of Fig. O are clear signatures 
of reverberation due to reflection. The width of the line-like 
feature in the lag spectrum does not necessarily mean that 
the actual spectral line is broad. A distant reflector with a 
narrow line at 6.4 keV and an edge at 7.1 keV can possibly 
prod uce a broad line-lik e structure similar to that seen here 
(e.g. iKotov et al.l l200ll for black hole binaries modelling) . 
Therefore, could this be due to the cold reflection compo- 
nent, which is responsible for the strong narrow line at 6.4 
keV in the flux spectrum (Fig. [T])? 

There are several points that can help answer the ques- 
tion. The change in the lag spectrum between the low and 
high flux states is dramatic. Fig. |31-top indicates that the lag 
is scaled by a factor of ^ 20, while the source flux changes 
by a factor of '-^ 4 — 5 (Fig. [l]). The flux of the narrow line 
in the flux spectrum does not change significantly with flux 
jSchurchet al.ll2003l ). so if the lags are produced in the nar- 
row line, the lag change with flux has to be due a change 
in the illuminating source and not in the reflecting medium 
or its geometry. In other words, a constant narrow line with 
flux cannot be responsible for a variable lag spectrum. 

Also, the value of the lag is too small to originate in 
the cold distant reflector. Early observations of the nar- 
row line variability constrained it to be at least a few light 
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Figure 5. The RMS spectrum calculated for the low and high 
flux observations. The difference resembles a reflection spectrum. 

years away from the cent ral source (e.g. iMarkowitz et al] 
I2OO3I : ISchurch et all |2003| ). Smaller lags can be produced 
by a distant refle ctor if the reflecting medium is patchy 
l|Miller et al.l l20ld'). However, the difficulty there is that the 
lag has an oscillatory behaviour that makes the lag measured 
over a wide frequency range, as is the case here, averages to 
zero. The alternative of using an energy-dependent trans- 
fer function could produce small lags over a wide frequency 
range but it has to be finely tuned by putting the refiecting 
medium exactly in o ur sight-line to produce the observations 
|ZoghbieiaL||20il). 

The lag spectrum in Fig. [3] has a wing that extends be- 
low 6.4 keV, which indicates a possible origin in a red-shifted 
line. Whichever component is producing the lag, it appears 
to be present (or more prominent) in the low flux spectrum 
and not in the high flux. Taking the difference spectrum be- 
tween the two states (e.g. observations 4 and 1 in Fig.©, 
we find that an absorbed power-law model is not sufficient, 
the data requires an additional component that matches the 
shape of the lag spec trum and resembles a broad iron line. 
The RMS spectrum l|Edelson et all I2OO2I ) for the low and 
high intervals is also different (Fig. [5}. At high fluxes, the 
RMS spectrum is decreasing smoothly with a drop at the 
6.4 keV, while there is extra variability between 4-6 keV at 
low fluxes. The shape of this extra variability component 
resembles a reflection spectrum. The RMS spectrum further 
confirms that there is a refiection component that is more 
prominent in the low flux observations. 

In order to further understand the differences between 
the low and high flux intervals, we fitted the spectra with 
an absorbed power-law plus cold reflection modelled with 
the reflection table reflionx (Ross & Fabian 2005). Moti- 
vated by the shape of the lag spectrum, we also included a 
relativistically-blurred reflection component. Again we used 
reflionx to model the reflection spectrum, this time, con- 
volved with relativistic kernel kdblur to account for the 
blurring due to special and general relativistic effects. Fig. [6] 
shows the best fit for the low and high flux intervals (plot- 
ted for observations 1 and 4 as an example). In the high 
flux interval, the spectrum consists of an absorbed power- 
law {Nh = 5 X lO^'^cm"^), a cold distant reflector (dot- 
ted line) and relativistically-blurred refiection (dashed). At 
low fluxes, the relative contribution of the blurred refiection 
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Figure 6. The spectra of NGC 4151 in the low and high states 
(observations 1 and 4 in Fig. [T] respectively), with the best fitting 
model. The continuous line is the total model. The cyan dot- 
dashed line is a power-law, the red dotted line is a cold reflector, 
and the green dashed line is a relativistically blurred reflection 
model. The residuals of the fit are plotted in the bottom. 



is higher. In other words, the main difference between the 
two flux observations is a difference in the reflection frac- 
tion from the blurred component. The reflection fraction 
is plotted in Fig. (Sjbottom, where it is taken as the ratio 
of the relativistically-blurred component to the power-law. 
This gives the relative contribution of each component to the 
different energy bands. The small residuals in Fig. [6] around 
6 keV are possibly due to the fact that we used reflionx to 
model the cold reflection. The model has a non-zero lower 
ionisation limit. 

Comparing the top and bottom of Fig. [31 the absence 
of lag in the high flux observations is easily understood. The 
observed lag has contributions from both the reverberation 
lag between the power-law and the reflection component as 
well as any lag between the power-law component and itself 
in the different energy bands. The strength of the contribu- 
tion to the observed lag from each component is determined 
by which component is dominating the variability. At high 
fluxes, the contribution of refiection is much smaller, and 
thus the measured lag will be dominated by the lag between 
the power-law in each band (equal ~ 0) rather than the 
reverberation lag. 

It is very interesting to note that the peak in the lag 
spectrum is at lower energies compared to the broad iron 
line (4.9 vs 6.4 keV). The two vertical lines in Fig. [3] mark 
the positions of the peaks in the lag and refiection fraction 
spectra. If only low frequencies are considered (Fig. [2]), the 
peak in the lag spectrum moves to higher energies. To ex- 
plore the apparent frequency-dependence of the lag spectra 
we extract lag spectra at two frequency bins < 2 x 10""" and 
(5 - 50) X 10"^ Hz (Fig. EJ. Then we fitted it with a sim- 
ple model consisting of a gaussian plus a constant. The two 
gaussians have energies 5.7 ± 0.19 and 4.60 ± 0.43 keV for 
the two frequency bands respectively. Confidence contours 
for the line energy and width are shown in Fig [7] Although 
there is an overlap in the width considering the 99 percent 
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Figure 7. Right: Lag-energy plots (similar to Fig. I2I3I I. Blue 
triangels are for < 2 X 10~^ Hz and red diamonds are for 
(5 — 50) X 10~^ Hz. The position of zero in this panel is arbi- 
trary, left Confidence contours for the width and energy of the 
gaussian lines fitted to the lag spectra. Dashed and continuous 
lines are for frequencies < 2 X 10~^ Hz and (5 — 50) X 10~^ Hz 
respectively. The contours represent 90, 95 and 99 percent limits 
shown in blue, red and green respectively. 

confidence contours, the energy is clearly higher at lower 
frequencies. This result has been tested also by calculating 
the lag spectrum at low frequencies, and accumulating more 
and more frequencies. As higher frequencies are added, the 
width of the line increases and the gaussian energy shifts 
to lower energies. This makes the result at low frequencies 
more robust and not just an artefact of the small number 
of averaged frequencies (see. end of sec.[2|. Interpreting this 
behaviour is straightforward if the delays are due to the re- 
verberating broad Fe-K line. Redder parts of the broad line 
(i.e. lower energies) originate in smaller radii and would nat- 
urally vary on smaller time-scales. Bluer parts of the line are 
emitted slightly further out and so vary on slightly longer 
time-scales. 



4 DISCUSSION AND CONCLUSIONS 

All previous lag detections in AGN have been in the soft 
band (< 1 keV). In this work we reported the first detec- 
tion of inter-band delays in the Fe-K band in an AGN. The 
lag is seen prominently in low flux observations and it is 
much smaller at high fluxes. The lag is in a sense that the 

5 — 6 keV band lags both 2 — 3 and 8 — 9 keV bands by 
~ 2000 seconds (this value is frequency-dependent). The 
interpretation of these delays is consistent, to first order, 
with a picture in which an illuminating power-law compo- 
nent, dominating at <~ 3 and > 7 keV is refiected by an 
accretion disk at distances of ~ 10 rg (The light crossing 
time at 1 rg is 5 x 10~^A/ seconds, where M is the black 
hole mass in solar units). The measured delays are due to 
light-travel time between the illuminating corona and a re- 
flecting accretion disk. The main difference between low and 
high flux observations is a difference in the reflection frac- 
tion. Because the energy bands used do not contain a single 
component (direct power-law only or reflection only), the 
reduced reflection fraction causes the delays to be smaller in 
the high flux observations because they are diluted by the 



direct component where there is no clear lag between the 
energy bands. 

The shape of the lag spectrum resembles a broad iron 
line peaking around 5 keV. Although the flux spectrum has 
a strong narrow line component, several arguments indicate 
that it is not the origin of the delays. First, the delays are 
small in value (~ 2000 s), which is much smaller than the 
distance scales of the region producing the narrow line (the 
molecular torus or the broad line region). Earlier work, con- 
strained th e region producing the narrow line to be on scales 
light years (|Ogle et al.l (|2000l )). The lag spectrum peaks well 
below 6.4 keV, or in other words, the most lagged band is 
not 6.4 keV corresponding to the peak of the line as would 
be expected if the narrow component is producing the lag. 
Furthermore, the shift of the peak with frequency (Fig. [7| 
rules out models where the del ays are due to path -difference 
close to the line-of-sight (e.g. iMiller et al.ll2O10l ). In those 
models, no frequency-dependent behaviour is expected, as 
the shape of the distant reflection spectrum is not expected 
to change on these time-scales. 

The light-crossing time-scale for NGC 4151 with a mass 
of 4 X 10^ Mq at say 6 rg is ~ 1200 s, which is of the same 
magnitude as those measured here. T he thermal time-scale 
at this radius for a standard thin disk (|Shakura fc SvunvaevI 
il973il . which could be the variability time-scale if ionisa- 
tion/density changes are involved, is ~ 2 x 10* seconds as- 
suming a viscosity parameter a of 0.1. Again, this time-scale 
is also consistent with the observations. Spectral fitting with 
a relativistic model gives an inner disk radius of 1.4 ± 0.1 
rg with an emissivity index of 4 ± 0.3, a profile flatter than 
usually observed. If the lag spectra are fltted with a rela- 
tivistic line model using the inclination and emissivity from 
the flux spectra, we find that the low and high frequency 
variability originate at ~ 10 and ~ 4 rg respectively. This 
also presumably indicate that the height of the illuminating 
source is not very large (a few rg ) . 

The lag spectrum of NGC 4151 consists of broad line- 
like feature that peaks between 4-6 keV. We have found 
that this peak shifts closer to 6.4 keV (i.e. higher energies) 
when lower frequencies are probed. This is consistent with 
the redder parts of the line emitted at smaller radii, and 
so they vary on smaller time-scales, and their light cross- 
ing time is smaller. Building an exact model however, re- 
quires a full treatment of the problem, including general 
relativistic effects and the interplay between the direct and 
reflected components in the observed data. To understand 
the lags, we plan to explore these details in the near future. 
This important result highlights what can be achieved with 
current telescopes in the brightest AGN NGC 4151. Larger 
collecting areas in future X-ray missions (e.g. LOFT and 
ATHENA), will make it possible to perform detailed studies 
in frequency, lag and energy parameter spaces in this and 
many other objects. 
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